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Abstract 

We compute the quasinormal frequencies corresponding to the scalar sector of gravi- 
tational perturbations in the four-dimensional AdS-Schwarzschild black hole by using the 
master field formalism of [1]. We argue that the non-deformation of the boundary metric 
favors a Robin boundary condition on the master field over the usual Dirichlet boundary 
condition mostly used in the literature. Using this Robin boundary condition we find a 
family of low-lying modes, whose frequencies match closely with predictions from linearized 
hydrodynamics on the boundary. In addition to the low-lying modes, we also see the usual 
sequence of modes with frequencies almost following an arithmetic progression. 
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1 Introduction 

Recently, a lot of attention has been devoted to the study of quasinormal perturbations 
in asymptotically anti-de Sitter (AdS) backgrounds. The first quasinormal mode (QNM) 
computation in AdS space was done in |2] for a conformally invariant scalar field, and then 
the problem was solved in for any minimally-coupled scalar field in dimensions d = 4, 5, 
and 7. The gravitational perturbations of global yl(iS'4-Schwarzschild, which is what we are 
interested in, have been computed for the first time in |3] . Since then, various properties and 
generalizations of these QNM's have been considered, such as asymptotic relations [H [6], 
difi^erent anti-de Sitter backgrounds [3 [8], or other boundary conditions [9]. 

We will phrase the problem in terms of the master field formalism that was developed 
by Kodama and Ishibashi in [T]. Based on previous ideas developed by Regge, Wheeler and 
ZeriUi [IHIE] that were further extended in [ll[7l[12], Kodama and Ishibashi developed this 
formalism to decouple the linearized Einstein equations for the gravitational perturbations 
of the global AdS-Schwarzschild in a gauge-invariant way and in any number of dimensions. 
In general, the perturbations in AdSa can be divided into tensor, vector, and scalar pertur- 
bations, depending on whether they correspond to expansions in tensor, vector, or scalar 
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spherical harmonics on the S**^"^ section of AdSa- The basic idea then is that we can ex- 
press each of these perturbations in terms of a master field $ and the appropriate spherical 
harmonics, and all we have to do is solve a certain differential equation satisfied by $. This 
equation will in general depend on both the perturbation type we are considering and on 
the number of dimensions. 

We will restrict our attention to the scalar sector of the perturbations in d = 4, where 
most QNM-related computations in the literature use a Dirichlet boundary condition on the 
master field $ near the boundary of AdS. The purpose of this paper is to comment on the 
choice of this boundary condition, and to suggest that a Robin boundary conditioiu would 
be more appropriate, especially from the point of view of the AdS/CFT duality ( [131 [HI [15] ! 
for a review, see [16]). It follows from the AdS/CFT dictionary that a natural expectation 
is to demand that the perturbations do not deform the metric on the boundary of AdS, 
and this condition in turn determines the asymptotic behavior of the master field at the 
boundary. While having no boundary deformations amounts in other similar situations to 
imposing a Dirichlet boundary condition on $ at the boundary, this is not the case for the 
scalar sector of gravitational perturbations in AdS^, where a Robin boundary condition is 
required (see section [3]) o Using the Robin boundary condition proposed in section [3l we 
find a family of low-lying modes that were not seen when a Dirichlet boundary condition 
was used instead (see for example [1]). In addition to the low-lying modes, we also find a 
tower of modes that is similar to the tower of modes found in [3] . For details on what makes 
the low-lying modes different, or for what other differences we find between our QNM's and 
the ones computed in [3], see section HI It is important to note that our Robin boundary 
condition doesn't affect the vector gravitational perturbations in AdS4^, because in this case 
the Dirichlet boundary condition is still appropriate. 

A good check on the values of the low-lying quasinormal frequencies comes from a lin- 
earized hydrodynamics approximation on S"^ x R. The rationale of this approach lies in 
the observation that since M-theory on y4(iS'4-SchwarzschildxS'^ is dual to a thermal CFT 
on the boundary, some QNM's should correspond to hydrodynamic modes of the thermal 
CFT. This idea has been developed in several interesting papers: in [T7l[T8l[T9] it was shown 



Robin boundary condition specifies a linear combination of a function and its derivative at the bound- 
ary. 

^We can anticipate some trouble in the scalar sector of AdS4 perturbations just by looking at the general 
large p dependence of the master field $ for any kind of perturbations and in any number of dimensions. 
In general, $ satisfies a second order differential equation whose linearly independent solutions behave like 
p~ and p~ ^ , where j ~ 0, 1, or 2 for scalar, vector, or tensor perturbations, respectively. The scalar 
perturbations in d = 4 are the first ones for which the behavior of the first family is subleading to the 
behavior of the second family. 
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that the quasinormal frequencies should correspond to poles of the correlation functions on 
the field theory side, and in [TSl El EO] this result was checked by explicit numerical compu- 
tations. We follow the approach in [21], where it was noted that the low- lying scalar and 
vector modes in /ive-dimensional AdS-Schwarzschild can be computed through a linearized 
hydrodynamics approximation. Extending the argument given in [21] to any dimension, 
we derive an approximate formula for the low-lying scalar and vector modes in AdSd- We 
find excellent agreement between the numerically found low-lying modes (using our Robin 
boundary condition) and the linearized hydrodynamics prediction in d = A. 

The paper is organized as follows: in section [2] we present an overview of the general 
setup of our calculation, in section [3] we comment on the choice of boundary conditions and 
derive the boundary asymptotics for the master field in section H] we show the results of 
our numerical computation of the quasinormal frequencies of the global y4(iS'4-Schwarzschild 
solution, and finally, in section [5] we compare our results to what one would expect from the 
analysis of linearized hydrodynamics of a conformal plasma on S"^ x R. 

2 Setup of the calculation 

In this section we briefly review the setup of our calculation. The global 74(i5'4-Schwarzschild 
black hole solution is given by 



where dfll is the standard metric on the unit S"^, 

dQl = -f^jdy'dy^ = dO^ + sin^ edcj)^ . (2) 
This metric is a solution to the Einstein equations that follow from the action 

S = ^ [ d'x^(R+^) . (3) 



2k^ J V L\ 

The horizon radius of the black hole solution ([T]) is then the positive root of the equation 

Po = Ph(i + §). (4) 
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For future reference, the mass, entropy, and Hawking temperature of this black hole solution 
are: 

We are interested in linear perturbations of the background metric ([I]), of the form gab + S gab, 
that satisfy the linearized Einstein equations following from ([3]). The boundary conditions 
satisfied by these perturbations will be discussed in section [31 

The linearized equations satisfied by the perturbations can be solved by separation of 
variables. We assume Sgab ~ e~*'^'^$(p)S'ab(^, 0), where the functions Sab depend only on the 
angular variables on S*^, and can be written in terms of the spherical harmonics Yim{0,(f)) 
and generalizations thereof. The exact equations describing the scalar, vector, and tensor 
perturbations of the d-dimensional AdS-Schwarzschild background can be found in [Ij. As 
discussed previously, we will only focus on the scalar perturbations in the case d = 4. Using 
the notation in [2T], we split the coordinates = (r, p, 6, (p) into = (r, p) and = {6, (p). 
We denote by Vi the covariant derivative with respect to the metric (|2j) on S*^, and by Da 
the covariant derivatives with respect to the two-dimensional metric 

d4 = -fdr' + W / = l-- + 7^- (6) 
The equations describing the scalar perturbations then read: 



SgalS = fal3 §(6*, (t>) Sgai = pfa Si (6*, (j)) 

6g,, = 2p2 [Hl{t, p) 7,, §(^, 0) + Ht{t, p) S,,(^, , 



(7) 



= -^diS = ^VidjS + l-fijS (8) 

ks ki. 2 ^ ' 

H = m + 3w w = — m = kl-2 (g) 

Xa = J— i fa + T-daHr 

ks \ ks 

Fat3 = fat3 + D^Xp + DpX^ (10) 

F = HL + \HT + -{d''p)Xa 
I p 

F\ = D^Fafs = 2dfsF (11) 
Fa(^ = ^ [d^Op (pH^) - ^gapD.d'^ (pH^)^ (12) 
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(^D^d^ - $ = (13) 

Vs(p) = + (2 + 3w) + 9mw^ + 9w^ (2/ + 3w - 2)1 . (14) 

where S denotes any of the spherical harmonics Yi^ on S^, and kg is the corresponding 
eigenvalue of the laplacian: 

{Vid' + kl)^ = o s{e,<p) = Yi^{e,<p) kl = i{i + i). (15) 

It is worth noting that the above master field formulation is gauge invariant. So equations 
(l71)-(fT^ don't determine the perturbations 6gab uniquely: there is an implicit freedom of 
choosing four of these functions through a gauge transformation of the form 6gab ^Qab + 
Vaffe + Vfot'a, where this time Va denotes the covariant derivative with respect to the full 
four-dimensional metric ([I]), and Va are arbitrary functions. A small discussion of our gauge 
choice and the residual gauge freedom is included in the Appendix. 



3 Boundary conditions 

3.1 Boundary conditions at p = oo 

The question of what boundary conditions one should impose on the master field $ at the 
boundary of AdS does not have a well-established answer: most of the previous authors 
have set $(cx3) = (see, for example, [221 123], O [6] ) , but other boundary conditions have also 
been used (see, for example, |9])o As we shall see below, the AdS/CFT dictionary relating 
perturbations and expectation values of operators in the dual field theory might help clarify 
this point. 

From the AdS / CFT perspective, there are two independent behaviors of the metric per- 
turbations 5gab at large p: 6gab ~ p^, which corresponds to a deformation of the boundary 
metric, and Sgab ~ 1/p, which corresponds to a non-zero VEV of the stress-energy tensor 
in the boundary theory. In defining the quasinormal frequencies it is sensible to require 
that the metric perturbations do not change the boundary metric, so they only produce a 
non-zero VEV of the stress-energy tensor (Tab) of the thermal plasma on the boundary. This 
prescription is equivalent to requiring the quasinormal frequencies to correspond exactly to 
the poles of the correlation functions in the strongly coupled dual CFT in the planar limit 



^It is not clear to us how the boundary condition that we find is related to the one proposed in [9]. 
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(see for example [H]). With this in mind, the significant challenge is to find the relation 
between the asymptotic behaviors of 5gab and $, which is what we'll now turn to. 
At large p, the master equation f|T3l) takes the form 



2^ fc| 18pg 1 

L2 (A;|-2)2 



far 



0, 



(16) 



where we have assumed e~^^'^ behavior and denoted u = Q/L. Being a second order dif- 
ferential equation, equation (fT6!) has two linearly independent solutions. Their asymptotic 
behaviors at large p are given by: 



$ far(p) = e-*^^/^ 
$far(p) = e-*^^/^ 



and 



P 



p3 



:i7) 



As noted earlier, the boundary condition that has been mostly used in the literature is 
ip^^^ = 0. As we shall see shortly, this condition is not consistent with the idea that Sgab ~ 1/p 
is the only behavior allowed. To argue this, we choose to work in axial gauge {Sgpa = 0), and 
we derive the boundary condition on $ required by 6gab ~ 1/p- While we include a detailed 
and more complete derivation in the Appendix, we now present the simplest way of arriving 
at the proposed boundary condition. 

Setting L = 1, we can plug ffTTI) into equation f[T^ and obtain, for the F^p component 
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Using axial gauge and, as discussed above, assuming 6gab ~ 1/p, we have: 



frp = Lp = L = 



Hr 



O 



P 



By using ( iTOl) we can compute 



P^ 



fr = 



O 



P" 



F = e" 



kg 



O 



(19) 



(20) 
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This means that axial gauge and 5gah ~ 1/P force F^-p to behave as By comparing this 

behavior to the general expectation fllSp . we conclude that the 1/ p term in fll8l) must vanish: 

^'" + ^ = 0. (21) 

Thus we obtain a Robin boundary condition, involving the master field and its derivative. 

3.2 Boundary conditions at the horizon 

In contrast to the large p boundary conditions whose derivation was somewhat subtle and 
tedious, the horizon boundary conditions are straightforward, being based on the requirement 
that classical horizons don't radiate. So in appropriate coordinates, the perturbations near 
the horizon should take the form of an infalling wave. To make this explicit, we define the 
standard "tortoise" coordinate by 

which puts the master equation (I13p into the form 

[-d'^ + dl-Vs{p)]<l> = 0. (23) 

Here, p —>■ pn corresponds to r^, —oo. Noticing that Vs{ph) = 0, we can immediately 
see that the near horizon behavior of the two linearly independent solutions to the master 
equation are 

$ncar(p) = f/e-^^(^+''* + ^e-if^C--.)/^' . (24) 
The infalling boundary condition then means setting = 0. 

4 Numerical solutions 
4.1 Change of variables 

In order to solve the master equation (1131) numerically, it is convenient to recast it in terms 
of a different field ip{y), defined by factoring out the near horizon behavior of the master 
field <l>(p): 

^ = e-^^i^+r,)/L^(y) y=l-^. (25) 

p ^ ^ 
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Setting L = 1, we can plug this ansatz into the master equation (fT3!) to obtain the differential 
equation satisfied by ip. We obtain 



[s{y)dl + t{y)dy + u{y)\ij{y) = Q^ 



(26) 



where 



s{y)=K{y){l-yYf 



t{y) = K{y) 



d 



(27) 



u{y) = -K{y)plVs 
K{y) = -[l + kl + ?>p]j-?>y{l + p]j 

y 

Here, K{y) has been chosen so that s{y), t{y), and u{y) are polynomial expressions in y that 
don't have any common factor and that don't vanish for any y between and 1. 

The remaining challenge before we proceed to solve the differential equation (126|) is to 
translate the Robin boundary condition for $ fl^T]) into a boundary condition for ip. This 
can be done by writing the first two terms in the series expansion of (125|) at large p: 



$(p) 



~ e 



-in{T+p,) 



(28) 



We get: 



^'(1) = — 

PH 



3po 



k1 2 



^(1). 



(29) 



Of course, the near horizon boundary condition ^ = translates into 



m = 1 , (30) 

and we can now turn to describing the numerical techniques that we use to solve the differ- 
ential equation fl2Bl) with the boundary conditions fl2^ and (13D]) . 



4.2 Numerical method and results 

Following the method used in [21] for the computation of quasinormal frequencies of the 
scalar modes, we integrate the differential equation fl26|) in three steps: 1) we develop a 
series expansion around y = and evaluate it at ?/ = = |; 2) we integrate the differential 
equation numerically by using Mathematica's NDSolve from y = yi io y = y f (to be given 



8 



« (ii) (iii) 







Figure 1: The zeroes of s{y) represented as crosses in the complex plane. The red cross at 
1/ = denotes a simple zero, while the green crosses denote double zeroes. We use a series 
expansion in region (i), numerical integration in region (ii), and another series expansion in 
region (iii). 



below); and 3) we match our numerical solution onto a series expansion around y = 1, which 
is computed using the boundary condition ( l29l) . In doing the matching, we compute the 
Wronskian between the numerical solution and the analytical approximation near y = 1. 
The Wronskian vanishes only when the two functions are linearly dependent, i.e. when 
satisfies the boundary condition fl2Ul) at y = 1. 

In developing the series expansions, we should keep in mind that the series solutions are 
guaranteed to converge only when s{y), seen as a function of the complex variable y, doesn't 
vanish. It is easy to obtain the zeroes of s{y) by writing 



siy) = y{y - viYiv - y2f{y - y2f 



(31) 



where 



1/1 = 1 



1/2 = 1 



kg 



3(1 + Pi,) 



ipH\/^ + 3p 



(32) 



2{l + pl) 2{l + pl) 

It follows that the series expansion around y = converges on the whole interval between 
and 1 (though the convergence close to y = 1 might be slow, because of the nearby zero 
of s{y)). Similarly, the series expansion around y = 1 has a radius of convergence r equal 



to the minimum of |yi — 1| 
shown that a good value for yf was yj 



Si) 



Ph 



P'„ 



(see figure [T]). Experience has 



1 - r/4. 

Using the method described above, we computed the lowest nine quasinormal frequencies 
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freq\/ 


2 


3 


4 


5 


6 


no 


2.156 - 0.285i 


3.361 -0.3542 


4.487- 0.333 i 


5.561 - 0.298 i 


6.608 - 0.266 2 




3.463 - 2.573 z 


4.461 - 2.443 z 


5.528 - 2.271 z 


6.577- 2.106 z 


7.610 - 1.963 z 


^2 


5.230 - 4.942 z 


6.023 -4.791 z 


6.964- 4.571 z 


7.935 - 4.340 z 


8.910 -4.126Z 


^^3 


7.096 - 7.308 i 


7.757- 7.165 i 


8.592 - 6.942 i 


9.484- 6.685 i 


10.40 - 6.432 i 


^4 


9.002 - 9.670 i 


9.572 - 9.540 i 


10.32 - 9.327i 


11.15 - 9.0642 


12.00 - 8.794 i 


^^5 


10.93 - 12.03 i 


11.43 - 11. 9H 


12.12 - 11. 7H 


12.89 - 11.45i 


13.69 - 11.18i 


n^ 


12.86 - 14.39 i 


13.32 - 14.28 i 


13.95 - 14.09 i 


14.68 - 13.842 


15.44- 13.56 i 


^7 


14.81 - 16.74i 


15.23 - 16.64i 


15.82 - 16.47i 


16.50 - 16.23 i 


17.23 - 15.95 i 


ns 


16.76 - 19.102 


17.14 - 19.00 2 


17.70 - 18.84i 


18.35 - 18.6H 


19.04- 18.342 



Table 1: Frequencies of scalar quasinormal modes for pn = 1 in units where L = 1. 



heq\pH 


0.2 


1 


10 


100 


Qo 


2.793 - 0.0008 z 


2.156 - 0.285 z 


1.739 - 0.066 z 


1.732 - 0.007i 




4.201 - 0.084 z 


3.463 - 2.573 z 


18.66 - 26.63 z 


185.0 - 266.4 z 




5.468 - 0.523 i 


5.230 - 4.942 i 


31.84 -49.172 


316.1-491.62 




6.896 - 1.121 z 


7.096 - 7.308 z 


44.95 - 71.70 z 


446.5 - 716.8 z 




8.416 - 1.735 i 


9.002 - 9.670 i 


58.03 -94.22 i 


576.6 - 941.8 2 




9.984- 2.348 i 


10.93 - 12.03 i 


71.10 - 116.72 


706.6 - 1167i 


fie 


11.58 - 2.956i 


12.86 - 14.39 i 


84.18 - 139.2 i 


836.6 - 1392 i 




13.20 - 3.559 i 


14.81 - 16.74i 


97.25 - 161.82 


966.5 - 1617i 




14.83 - 4.158 i 


16.76 - 19.10i 


110.3-184.32 


1096 - 1842 i 



Table 2: Frequencies of scalar quasinormal modes for / = 2 in units where L = 1. 

for ph = 1 and / = 2, 3, 4, 5, and 6 (see table [1]), and for pn = 0.2, 1, 10, and 100 at fixed 
/ = 2 (see table [2]). In these tables we have included only the quasinormal modes with 
Refi > 0; equation 0261) implies that if fi is a quasinormal frequency, then so is — fi*, so the 
QNM's with negative real parts can be trivially obtained from the ones with positive Refi. 

The most prominent feature of the quasinormal modes included in tables [T] and [2] is the 
separation of the quasinormal frequencies Qn into two groups: a main series of fast modes 
given by Qn with n > 1, and low- lying slow modes given by Qq (for a similar feature of 
the quasinormal frequencies in ArfS's-Schwarzschild see [21]). The low- lying modes differ 
significantly from the fast ones in a number of ways: 

• While the fast modes form a tower of modes at each value of ph and /, the low-lying 
modes stand out as not being part of this tower (see figure [2]). 

• The low-lying modes have a different p/^-scaling from the main-series ones (see figure[2]). 



10 





5 


10 


20 


50 


100 


2 


1.761 - 0.129 z 


1.739 - 0.066 z 


1.734 - 0.033 z 


1.732 - 0.013 z 


1.732 - 0.007 z 


3 


2.547- 0.312 z 


2.474- 0.164 z 


2.456 - 0.083 z 


2.450 - 0.033 z 


2.450 - 0.017z 


4 


3.380 - 0.537Z 


3.219 - 0.292 z 


3.177- 0.149 z 


3.165 - 0.060 z 


3.163 - 0.030 z 


5 


4.273 - 0. 787 z 


3.979 - 0.449 z 


3.900 - 0.231 z 


3.877- 0.093 z 


3.874- 0.047 z 


6 


5.230 - 1.043 z 


4.761 - 0.631 z 


4.628 - 0.329 z 


4.590 - 0.133 z 


4.584- 0.067 z 


7 


6.246 - 1.286 i 


5.566 - 0.836 i 


5.362 - 0.442 i 


5.303 - 0.180 i 


5.294- 0.090 i 


8 


7.311 - 1.499i 


6.399 - 1.059 i 


6.103 - 0.570 i 


6.017- 0.233i 


6.004 - 0.117i 


9 


8.410 - 1.675 i 


7.262 - 1.297i 


6.852 - 0.713 i 


6.731 - 0.292 i 


6.714 -0.147i 



Table 3: Frequencies of some of the low-lying modes in units where L = 1. 
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(b) PH = 1 



Im{Q} 



«» 


-100 -50 


50 100 


m-25 
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• -50 


• 


• -75 
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• -100 
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• 
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• 


-175 
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(C) PH - 10 



Figure 2: Quasinormal frequencies for pn = 0.2, pn = 1, and pn = 10, in units where L = 1. 
The black dots represent the main-series modes, while the red ones represent the low-lying 
modes. It is fairly clear that for pn = 0.2 and 1 the low-lying modes are not part of main 
series tower. This is not obvious in the pn = 10 case, because of the plot scale. 
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Figure 3: Quasinormal frequencies for different values of ph plotted against /, in units where 
L = 1. The blue stars correspond to pn = 5, the red triangles to pn = 10, the dark green 
diamonds to pn = 20, the light green triangles (barely visible) to pn = 50, and the dark 
blue stars to pn = 100. The dotted line represents the linearized hydrodynamics prediction 
f l33|l . which matches almost perfectly the numerical results for large pn- 



This feature is most clearly seen at large pn, where the low- lying modes approach 
r2 = a//(/ + l)/"\/2 as ^ oo (see next point), while the main series modes grow 
proportional to pn'- compare, for example, the columns corresponding to pn = 10 and 
Ph = 100 in tabled 

The low-lying modes can be interpreted as the linearized hydrodynamic modes of a 
conformal plasma on S*^ x R. While we will explain this correspondence in more 
detail in section \5[ for now it is worth mentioning that a linearized hydrodynamics 
approximation on S'^ x R gives, up to first order in 1/ph, that 



GPH 



Ph 



(33) 



with ks = + 1). A plot of low-lying modes for various values of / and pui together 
with the hydrodynamics prediction ([3SD can be seen in figure [31 which is based on the 
numerical values in table [31 

It is worth noting that while the tower-like feature of the scalar QNM's can be observed 
even if one imposes a Dirichlet boundary condition, the low-lying modes have not been 
seen in either numerical computations or analytical approximations that use the Dirichlet 
boundary condition (see, for example, [221 [6]). 

Leaving the low-lying modes aside, we can compare the structure of the main series 
modes to the structure of the modes described in [22] that come from imposing the Dirichlet 
boundary condition on the master field. We find that the spacing between the main series 
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modes at large f2 asymptotically approaches the spacing between the modes computed in 
[1]. However, the initial offset of the tower is different, our modes being in between the ones 
found in [1]. 



5 Linearized hydrodynamics 

In [21] it was noticed that in the case of the global AdS^ black hole there was a separation in 
the imaginary parts of low-lying scalar modes and the "main series" modes. The former were 
interpreted as hydrodynamic modes and the latter as microscopic. So a simple treatment 
of linearized hydrodynamics should be able to reproduce these low-lying modes in other 
dimensions as well. The goal of this section is to develop such a treatment. 

In thinking about hydrodynamics, the general setup on S'^~^ x R is given by the following 
relations: 

Tab = (e + P)UaUb + PQab + Tab 



Tab = -V [ AacVUb + AfocV^a " j^AabV^U^ ) - ^A^iVSc 

Aab = gab + UaUb 
y^Tab = . 



(34) 



where cjab is the metric on S'^~'^ x R and Vq is the covariant derivative with respect to this 
metric. Since the theory on the boundary of AdS is conformal we expect T\ = 0, which 
implies both e = {d — 2)p and = 0. Following the same approach as in |2Tj, we ignore 
the temperature- independent contribution from the Casimir energy to T"''. The Casimir 
energy comes from considering the quantum field theory on the compact space S'"'"^. For 
our purposes we can think of it as a temperature-independent shift of the zero point energy, 
which can be safely ignored. 

The vector describes the velocity at each point in the fluid, and we choose to normalize 
it by imposing u'^Ua = —1. Let us denote u"" = (1,m*) where i runs over the S'^"^ directions. 
In the linearized approximation we consider to be small. Perturbing at the same time the 
pressure p = po + Sp, one can derive from flM|) the linearized equations 



(35) 



{d-2)^ + id-l)poV,u' = 



{d - l)Po^ + V*(5p - 7] (v^Vju' + VV^ j - Vj^d'VjU^ = . 
Note that for d = 5 equation (135|1 reduces to the linearized Navier-Stokes equation on 
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which was analyzed in section 5.3 of [21]. We wish to examine scalar perturbations next, 
which are described by the ansatz 



Sp = Kie-'^^E u' = i^2e-'^"V'§ , (36) 

where S satisfies ^V^ + S = 0, as explained in section O and L = 1. Plugging (l36l) into 
(153]) . we obtain the following system of equations for Ki and K2: 

-in{d - 2)Ki -{d- l)ksPoK2 = 

/ rf-3 \ (37) 

-ksK, + l-tn{d - l)po + '^Vjzr^kl - 2{d - 3)r]\ K2 = . 

In order to have non-trivial solutions, this system must have zero determinant. This gives 
a quadratic equation for Q, whose solutions can be given in terms of a series expansion in 
V/Po- 

n_, ks ^v kUd-3)-id-2){d-3) fv'\ 

~ ~ Vo {d-i){d-2) + ^ U,) 

We can connect this result to the AdS^, quasinormal mode problem by noting that 



Po s 1 + ' 

which can be easily derived from ([5]) in the case c? = 4, but it is actually true in any number 
of dimensions. Using the conjectured lower bound on the viscosity 7 = ^ [211 [25], that has 
been checked in the AdS^ case in [26], we find 

v/d32 p„ (d-\)(d-2) \pIJ * ' 

It is easily seen that this reproduces the hydrodynamical modes previously discussed in the 
global AdS^ black hole case in [21]. For d = 4, equation (140|) reduces to (1331) . 
Similarly, we can describe the low-lying vector modes by the ansatz: 

Sp = u' = K^e-'^^V' . (41) 
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We find that the corresponding frequencies f2 are given by 



i{d-l) 




(42) 



It is interesting to note that the numerical value given by this formula when d = 4, I = 2 
and ph = 100 agrees within 10% with the low-lying vector mode of Table 9 in [22] . 

6 Conclusions 

In this note we examine the relation between the asymptotic behavior of the master field 
and the behavior of the scalar sector of metric perturbations in the global AdS4^ black hole. 
We argue that the boundary condition that corresponds to a non-deformation of the metric 
on the boundary translates into a Robin boundary condition for the master field. With 
this boundary condition, we compute the scalar quasinormal modes. We find some low- 
lying modes that have not been previously observed, and compare them with the linearized 
hydrodynamical modes of the boundary CFT. 
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A Gauge freedom and boundary conditions 



In this section we derive the asymptotic expressions for the metric perturbations 6gab by 
solving the system of equations (171)- (fT3l) . In doing so, it is important to reahze that equations 
(l711-( fT3l) don't determine the metric perturbations uniquely. We have seen that the gauge 
freedoms 

Sgab Sgab + "VaVb + V bVa (43) 

present in any perturbation theory problem in general relativity enables us to set 5gpa = 
(this is what we referred to as axial gauge). However, even after we set Sgap = we still have 
a residual gauge freedom left, and we would like to understand this residual gauge freedom 
a bit better before we derive the asymptotic expressions for 6gab- 

The first thing to note is that generic gauge transformations of type fH3|) do not preserve 
the form of the metric ([7]). Instead, generic transformations would just map our initial 
solution onto something that doesn't transform under the 5*0(3) isometry group of S"^ in 
any definite way. We only look at perturbations with the specific 5*0(3) structure defined 
in ([7]). Therefore, we need to restrict the class of allowed gauge transformations to the ones 
that preserve this 5*0(3) structure. Such transformations are of the form 



VaiT,p,e,(l))= (K{t,pP{9,^) hp{T,p)S{9,^) h{T,p)MOA) h{T,p)S^{9,^)) , (44) 
and they give 



2 V(it;j) = -2hks^,, + {hks + 2pfhp) 7i,§ . 



dc,h - haks - -{dap)h 
P 



(45) 



It is easy to see now that if we start with any scalar perturbation of the form (JTj), we can set 
^Opa = by solving three first order non- homogeneous differential equations for hr, hp, and 
h. The residual freedom that remains after setting Sgpa = is reflected in the choice of the 
three integration constants (which are functions of r) that enter in the general solutions of 
these equations. So in addition to setting 6gpa = we also have the freedom to prescribe the 
time behavior of three of the other components of Sgab at a given point. In particular, the 
gauge freedom allows us to set the large p behavior of three such components to have no p"^ 
term in a large p series expansion. Requiring all of these components (which are described by 



"^In this section Va denotes the covariant derivative with respect to the four-dimcnsional metric ([T]). 
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the four functions Ht, Hl, f-^r) to have no terms cannot be accomphshed by making 
a gauge transformation, and is therefore meaningful as a boundary condition on the metric 
perturbations. 

We now turn to the problem of finding the asymptotic expressions for the metric coef- 
ficients 6gab in axial gauge. We will set L = 1 throughout the entire calculation. With the 
assumption 



fpp fp 







fr 



-iflr 



B{p) fr = e-^''^C{p) 



we first compute the quantities F and Fap that enter in equations ffTTl) and f|T2|) : 

^ ^''\kl[2AL{p) + AT{p)]+2pfA'M) 



2kg 



TP 



pp 



{-2zkspnC{p) + klB{p) - p2 [2n'AT{p) + ff'A'M]) 
(ksfCip) - 2zQpf [At{p) + pA'M] + 

+ pfksC'ip) + zfipV'^r(p) - kspf'Cip)) 
(A^(p)[4/ + p/']+2p/A^(p)) , 



(46) 



(47) 



and FpT- = F^-p. The plan now is to plug the above expressions into equations f|TT]) and f|T2|) . 
and find a series solution for the corresponding differential equations. In order to do this 
though, we need to get hold of the right-hand side of equation (fT2!) . perhaps in the form of 
a large p series expansion. This can be done by solving the master equation (fT3l) : 



where 



$(p) = if 



(1) 



P 



(2) 



(3) 



P' 



+ 



(48) 



(2) 



(3) 



9pi 



fi2 



(fc| - 2)2 
18pg 



(0) 



Po(2 + A:|) 

(fc| - 2)3 2{kl - 2) 



3pg 



(A;5 - 2) 



6 



(1) 



(49) 



and all higher order terms can be expressed in terms of linear combinations of y^*^''-* and ip'^^\ 
The two constants tp^^'^ and tp^^^ can thus be interpreted as the two integration constants that 
appear when we integrate the master equation, which is a second order ODE. The above 



17 



expansion can then be used to find a series expansion of the right-hand side of equation (fT2|) . 
The resulting expressions are long and not that insightful, so we will not reproduce them 
here; their derivation is nevertheless straightforward. 

We now solve for the functions A^^p), At{p), B{p), and C{p) that completely determine 
the metric perturbations via 



pi 







(50) 



in four steps: 

1. We first solve for At{p) from the Fpp equation in f|T2|) with the LHS given by the 
corresponding expression in (j47j) and the RHS computed from the series expansion 
We find: 



At{p) = A 



(0) 



A. 



(2) 



+ 



ksPo 



2{kl - 2) 



(1) 



kl-2) 



+ 



^ s 

12 



+ 



(0) 



1 / 1 



(51) 



Here, we can think of A^^ and A'^^ as integration constants: we have two integration 
constants because the differential equation satisfied by At{p) is second order, as can 
be easily seen from the Fpp relation in (147|) . 

2. Next, we solve for C{p) from the equation in fll2p . Again, the LHS of this equation 
is given in fHTl) . and the RHS can be computed from (jUj). We obtain: 



C(p)=C(-i)p + 



Spoip 



(0) 



/c| 2 



ksC^-''^ +iQ{2AP - aP] 



ks 



1 
P 



PoitnA^ - ksC(~^^) iksn{kl-2) 

' 



ks 



6 



4 + 

P VP 



(52) 



where C*-"^-' is an integration constant — we can see from the Frp relation in ( 1471) that 
the corresponding differential equation for C{p) is a first order ODE, so its solution 
has to have one integration constant. 



3. Similarly, we next solve for B[p) from the Frr equation in (lT2l) . As can be seen from 
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the Ft-t- relation in (H7|) . this equation doesn't involve any derivatives of B{p), so its 
solution doesn't involve additional integration constants: 



kl-2 + 2n^ 



+ 



2Q{iksC^~^'^ + QAP) - 2(1 + fi2)A 



A{ks 
+ 



2 -21 



(1) 



(0) 



kl 2 



, (1) , ^Pof 



(0) 



2pon{iksC'-~^'^ + 1]A 



(0)^ 



^5(As__^ (0) 

6 ^ 



1 / 1 

- + ^ - 

P VP 



(53) 



4. Finally, we solve for Al{p) from the second equation in ( fTTl) with (3 = t. Again, this 
equation doesn't involve any derivatives, so we have no integration constants: 



Al{p) 



4A 



(2) 



2kg 



+ 



+ 



{kg 2)^^-*^ ^ kl 2 



/c| 2 / 



p2 



2A;| 



(54) 



/c| 2 / 
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It can be checked that the above series solutions automatically satisfy the other two 
equations in fill I) that we have not used. Also, the fact that the integration constants A^j^\ 
Af, and are still undetermined is a consequence of the residual gauge freedom that 

we discussed at the beginning of this section: these three integration constants allow us to 
set the values of three of the functions Al(p), At{p), B{p), and C(p) at a given point to 
whatever we want. 

Requiring that the metric perturbations don't grow like at large p and using (!50|) - (!54|) . 
we get aP = A^ = C^"^) = 0, together with 



(1) I 3poy(°) 
kg — 2 



0, 



(55) 



which is the same as (ET 



We note that the relations ( l55ll and A^ = A^ = C*^ = make almost all terms 
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written in the series expansions (l5T|) - (l54l) disappear, and we're left just with 

_ kUki-2) y>' (i\ 



p \p^ 



.^MMn^e-^-^ + o(- 
6 



which looks incredibly similar to the expressions found in section 3.3.2 of [21] in AdS^- 
Schwarzschild. In light of the analysis done in [21], it is worth mentioning that the leading 
coefficients in ( |56|) give, up to a proportionality factor, the expectation value of the stress- 
energy tensor in the boundary 2 + 1-dimensional CFT. Conservation and tracelessness of the 
stress-energy tensor can then be easily checked using the same approach as in [2T] . 
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